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Abstract 

An interesting analog circuit for simulating a signal with fluctuations having 
a probability density function with a power tail has recently been proposed 
and constructed. The exponent of the power law can be fixed by tuning an 
appropriate circuit element. The proposal is to use the circuit as a simulator- 
generator of financial market fluctuations and as a tool for risk estimations 
and forecasts. We present a discussion of the stability conditions for multi- 
plicative noise and an exhaustive analysis of the power law fluctuation gen- 
erator in connection with the electronic components and their parameters. 
From our studies one can conclude that the proposal is not adequate to 
provide a confident experimental scheme to follow the fluctuations of the fi- 
nancial market due to both electronic implementation and to difficulties in 
parameter tuning. 
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1 Introduction 



Recently, Sato, Takayasu and Sawada have constructed a very interesting 
analog circuit able to generate a signal corresponding to a time series with 
fluctuations resembling those of a typical financial market. The correspond- 
ing probability density function presents a power law tail. 

It is worth mentioning that the proposal of Ref. |J immediately attracted 
the opinion of the scientific community. In particular, we should mention 
that the American Institute of Physics Bulletin of Physics News number 
478 of April 6, 2000 includes a comment on Ref. [|1] under the title "An 
electrical circuit mimics yen-dollar fluctuations". It ends by saying that the 
circuit costing approximately $5 can estimate yen-dollar fluctuations as fast 
as the $10, 000 workstations that are running mathematical simulations of 
the exchange rates. On the other hand, CERN Courier of June 2000 includes 
also a comment under the title "How can you make a million?" where one 
can read that Japanese Scientists have unveiled a $5 electrical circuit that 
can mimic fluctuations in the yen-dollar exchange rate. 

The circuit constructed in Ref. |l| being very interesting in its own, 
deserves a careful analysis from the electronic point of view. We present 
here an exhaustive study of the circuit behavior and of its potentialities for 
reproducing financial market fluctuations. We conclude that the circuit really 
mimics these fluctuations but the meaning of this mimics is precisely the one 
the Oxford Dictionary gives: imitate closely in order to make fun of. 

Before the presentation of the circuit and its relationship with the Langevin 
equation with multiplicative noise, we discuss, in Section 2 the stability con- 
ditions for this kind of noise. Afterwards we include in Section 3 an analysis 
of the power law fluctuations generator in connections with the electronic 
components and their corresponding parameters. A circuit simulation com- 
pletes our study and its results are presented in Section 4. Finally, Section 5 
summarizes our conclusions. 



Before entering into the discussion of both the electronic and the simulation 
analysis, let us briefly review the stability in multiplicative stochastic pro- 
cesses. It is well known that an additive stochastic process is stable whenever 
the associated deterministic problem is. In fact, given the process defined by 



2 Stability 
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where r] is a standard Gaussian noise, it has a stable stationary behaviour 
when all its moments (v n ) exist up to the n-th order. This imposes a condition 
on the integral of L(v) @. 

However, for a multiplicative noise the stability of the deterministic prob- 
lem is not enough to ensure the stability under fluctuations |2| . Consider now 
the process 

dv , , s 

— = —dv + vri (2) 
dt 

where rj is again a standard Gaussian noise and d a parameter. In this case 
the moments are given by 

( v n ) = (v n ) t=0 exp [-n t {d-n D/2)] (3) 

where D is the measure of the Gaussian fluctuation, namely 

( V (t + r)r ] (t)) = D5(r) 

The stationary distribution results 

p (v) = v- 1 - 2d / D (4) 

which is not normalizable. It is then clear from Eq. @ that in this case of 
multiplicative noise, the condition 

2d 

has to be fulfilled in order to guarantee the finiteness of the moments (v n ) 
when t — > oo. Consequently, if d < D/2, the multiplicative noise overcome 
the deterministic restoring term —dv, even if this deterministic problem has 
a stable solution. 

Having these mathematical conditions in mind, we re-analyzed the circuit 
proposed in Ref PJ. The main point is to simulate, electronically or numeri- 
cally, the equation of Langevin type with multiplicative noise (0), that in the 
present case reads 

^ = h4 + ^ (t) ) W(t) (5) 

Notice a change of sign in the first term of the r.h.s of this equation, with 
respect to Eq. (12) of Ref.Q. 
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Figure 1: Block diagram of analog simulation 

3 Electronics 
3.1 Generalities 

The general block diagram in Fig. [I] represents the integrated form of the 
original differential equation |5|. 

The different blocks involved in Fig. [I] are implemented by using elec- 
tronic components, mainly amplifiers. From the diagram in Fig. [I], t ne 
electronic implementation can be easily derived giving rise to the circuit in 
Fig. | 
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Figure 2: Power fluctuation generator 
The electronic components have their own characteristics and impose 
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some restrictions on the functions to be implemented. It means that the range 
of values where the function is valid is bounded by the electronic components 
behaviour. 

Let us then analyze, following Ref. Q, the integrator block implementa- 
tion shown in Fig. [3[ 




, |A| [db] 
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A(ito) = Vo/Vi 




Figure 3: Electronic integrator and its Bode plot 

Considering an ideal amplifier, namely one with Gain Product Bandwith 
(GBW)^ oo and Gain (A ) — > oo, the corresponding transfer function is 

V (s) _ R f 1 



Vi(s) R v (1 + sR f Cf) 
with s = iu. Fig. [5] includes also the Bode plot of this expression. 
Under the condition uj ^> l/(C/i?/), Eq (|6]) gives rise to 



(6) 



and 



V n (s) = - 



v (t) 



Rv (sRfC f ) 
1 



Vds) 



Vi(t) dt 



(7) 



(8) 



R f c f 

Consequently the above circuit works as an integrator for frequencies above 
uj p = l/(Rf Cf). In summary, to implement an integrator using a real oper- 
ational amplifier (OA) its GBW must be larger than the highest frequency 
of interest and A 3> Rf/R v . 

To verify the above conditions, the Bode plot of the desired integrator, 
should be compared with that of the open loop gain of the OA. As a rule of 
thumb, the implementation will be valid if the plot of the OA "contains" the 
integrator's Bode plot. 
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3.2 Circuit analysis 

The criteria above was applied to the original circuit proposed in Ref. [Q 
that includes the values 

C f = 10 pF ; R f = 100 KVt ; R V = 5-200Q 

giving rise to 

Rf 

uj p = 1 MHz ; A = -i- = 54 - 86 db. 

R v 

In Fig. ^, a Bode diagram of Eq. (^) is plotted together with the OA 
LF157 characteristic obtained from the National Semiconductor Co. Linear 
Databook. As can easily be seen, the condition that the OA "contains" 
the desired characteristic is not fulfilled. So the changes of R v between the 
specified values do not produce a full change in the "integrator gain" and 
the observed effects are only marginal ones. In other words, a change in R v 
is not fully reflected as a change of the " integrator gain" . 




Figure 4: Comparison between Bode plots for the amplifier and the integrator 

Moreover, when the circuit in Fig. || was tested with the specified values 
of R v , the output presents unwanted saturation as shown in Fig. ||. 

Only for values of R v > 1.7 KQ, the amplifier remains almost linear 
giving rise to the output seen in Fig. |^. In fact, the reported behaviour in 



6 




Figure 5: Output of the circuit with the original parameters for R v = 35 f2 
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Figure 6: Output of the circuit with the original parameters for R v = 2200 f2 



Ref. |TJ was qualitative found for 1.7 KQ < R v < AKQ. For larger values, 
the electrical noise dominates. 

It is important to remark that our measurements were done with a Tek 
scope (Model TDS 3030) with 8 bits of resolution. Nevertheless we are certain 
that this technical detail doesn't change our conclusions. 

Following the circuit criteria discussed above, a new set of values for R v , 
Cf and Rf was chosen . To include the integrator characteristic inside the 
Operational Amplifier (OA) open loop, a lower pole frequency (u p ) was in 
order. To avoid output saturation, higher values of R v are also needed. Then, 
we decided to use the values 

C f = 11 pF ; R f = 800 KQ ; R v = 6 - 12KQ. 

Resulting in 

u p = 113 KHz ; A = = 35 - 45 db 

R v 

This is not an arbitrary selection because a value of Cf higher than 15 pF put 
the system in a heavy oscillation. A careful study of the LF157 characteristics 
shows that the stability conditions are greatly affected if a capacitor larger 
than 15 pF is connected between the input and the output. Consequently, a 
different amplifier should be used instead. 

The circuit with modified parameters worked as predicted. In Fig. [7| we 
present the log-log plot of the probability density function for the R v values 
specified. 

Using the adequate values for CfRf, a tighter control of the integrator 
gain is obtained. Notice also that for lower values of R v , that means higher 
gain, the circuit does not met the integrator condition, and for R v higher 
than 12 KQ the electrical noise dominates. 

4 Simulation 

In order to check the influence of the amplifier's offset voltage we have per- 
formed additional tests. This offset, independently of the circuit values, does 
not affect the final power law behaviour and can be reduced in the data 
processing stage. 

If the circuit is to be really useful in predicting fluctuations on a quan- 
titative way, it must be necessary to know the influence of each electrical 
parameter on the output signal in order to reproduce, for example, the prob- 
ability distribution of historical data. This work would be tedious, when not 
impossible using the real circuit. 
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Figure 7: Probability density function for the modified parameters 

We propose to use instead some simulation procedures that, without loos- 
ing the electrical point of view, put in evidence the influence of each of the 
parameters. We have tried several configurations using the Spice program [f|]. 
Although this program allows one to use realistic components, by means of 
transistor- level models, we prefer to use quite ideal representation where the 
different parameters of the circuit, desirable or not, can be added separately. 
It can be mentioned that the simulation with real components simulation 
is also plagued of serious convergence problems, being certainly difficult to 
make it runs in a fiable way. 

Fig. [8] shows the circuit used in the simulation. It reproduces the dif- 
ferential equation of interest. Different gain blocks allows one to change 
amplitudes, parameters such as offset voltage and saturation limits when 
added, as well as different noise generators. 

The typical output of the simulation is summarized in Fig. ||. 



5 Conclusions 

We have carefully analyzed the recently proposed electronic analog circuit 
that generates fluctuations with a probability density function having power 
law tails. 
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Figure 8: Simulated circuit 
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Figure 9: Output of the simulation. a)signal; b) density of probability 
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We found that the originally proposed parameters are not the most per- 
tinent ones because they place the circuit in the borderline of the integrator 
Bode plot. With appropriate parameters, the circuit performance is note- 
worthily improved. Moreover, by changing the integrated circuit originally 
chosen, a further improvement can be achieved. 

On the other hand, the numerical simulation, even if carried out in terms 
of ideal electronic components, reinforces the previous conclusions. 

Finally, several comments are in order: 

• Even if the power fluctuation generator circuit is certainly a low cost 
apparatus, a costly peripheral equipment is needed for taking and an- 
alyzing the output data. 

• The circuit, in order to mimic reasonable fluctuations, needs a fine 
tuning of the pertinent parameters. 

• A direct numerical simulation of the original equation seems to be a 
much better solution. 
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